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ABSTRACT 

A simple model of chemical enrichment in cluster early-type galaxies is presented where 
the main parameters driving the formation of the stellar component are reduced to 
four: infall timescale (tj), formation epoch (zp), star formation efficiency (Ces) and 
fraction of gas ejected in outflows (iJout)- We find that only variations in J5out or CeS 
can account for the colour-magnitude relation, so that the most luminous galaxies 
had low values of ejected gas and high efficiencies. Less massive galaxies can be either 
related to a lower star formation efficiency (Ccff-sequence) or to an increased outflow 
rate (Bout-sequence). The combination of chemical enrichment tracks with population 
synthesis models (Bruzual & Chariot 2000) is used to explore the correlation between 
mass-to-light ratios and masses. A significant slope mismatch is found between stellar 
and total M/L ratios, which cannot be explained by an age spread and implies a non- 
linear correlation between total and stellar mass: Mtot oc Mg^. The sequences driven 
by star formation efficiency (Ces) and outflows (Bout) are shown to predict different 
trends at high redshift. The variation with redshift of the slope of the fundamental 
plane will increase signiflcantly in the efficiency sequence — driven by age — and will 
slightly decrease in the outflow sequence — driven by metallicity. The evolution of 
the zero point is similar in both cases and within the observational errors of current 
observations. Measurement of the dependence of the tilt of the fundamental plane 
on redshift will break the degeneracy between outflows and star formation efficiency, 
which will enable us to determine whether the colour-magnitude relation is controlled 
by age or metallicity. 
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1 INTRODUCTION 

The process of galaxy formation and evolution can be ex- 
plored in two complementary ways: A "forwards" approach 
takes into account the physics underlying the most basic pro- 
cesses of structure evolution and star formation and — after 
finding a suitable set of initial conditions — evolves the sys- 
tem forward so that the final output is compared with obser- 
vations. This is the philosophy behind N-body simulations 
or semi-analytic modelling (e.g. Baugh et al. 1998; Kauff- 
mann & Chariot 1998). On the other hand, a "backwards" 
approach simplifies the physics behind galaxy formation and 
evolution to a phenomenological problem comprising a re- 
duced set of parameters, using local observations as con- 
straints. Then the system is evolved backwards so that the 
predictions for a given set of parameters are compared with 
observations at moderate-to-high redshifts (e.g. Bouwens, 
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Broadhurst & Silk 1998a,b). This "brute force" method ab 
lows one to search a reasonable volume of parameter space, 
throwing light on processes, such as star formation, which 
are otherwise too complicated to tackle from basic princi- 
ples. 

Cluster early-type galaxies are ideal candidates for a 
comparison between models and observations. Significant 
samples of these galaxies can be found over a large redshift 
range. Furthermore, the tight observed correlations such as 
the colour-magnitude relation or the fundamental plane can 
be used as powerful constraints on a phenomenological back- 
wards approach. The list of observed clusters at moderate to 
high redshift is quite large (e.g. Dressier et al. 1999 ; Stan- 
ford, Eisenhardt & Dickinson 1998; Oke, Postman & Lubin 
1998; Van Dokkum 1999) and ever increasing (e.g. Yee et al. 
1999). Clusters observed at redshifts z ~ 1 — 1.5 yield valu- 
able information about the epoch of star formation, pushing 
it to very high redshift {zf ^ 3). 

Unfortunately, the direct spectrophotometric determi- 
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nation of the star formation history is hampered by the age- 
metallicity degeneracy (Worthey 1994; Ferreras, Chariot & 
Silk 1999) which allows the variations of most of the spec- 
trophotometric observables to be explained either by a range 
of ages or metallicities. Broadband photometry is strongly 
dependent both on age and metallicity but even spectral 
indices targeting single lines such as Balmer absorption or 
magnesium abundance can change both with age and metal- 
licity. Not surprisingly, a similar age estimation technique 
based on spectral indices applied to similar sets of elliptical 
galaxies yields contradictory results: Kuntschner (2000) and 
Kuntscher & Davies (1998) find coeval stellar populations 
in Fornax cluster ellipticals so that the colour range is ex- 
plained by a metallicity sequence. On the other hand, the 
sample of field and group ellipticals observed by Gonzalez 
(1993) and further analysed by Trager et al. (2000) presents 
a relatively large spread in ages. Hence, the issue of the stel- 
lar age distribution in galaxies still requires the aid of mod- 
elling. We will show that incorporating chemical enrichment 
allows one to potentially solve the age-metallicity degener- 
acy problem. The next three sections describe our chemical 
enrichment model and the meaning of the reduced set of pa- 
rameters used to trace the star formation history in cluster 
ellipticals. §5 deals with the comparison of predicted and 
observed mass-to-light ratios and its use at high redshift 
to discriminate between a mass sequence driven by age or 
metallicity. Finally in §6 we discuss the predictions and list 
the main conclusions. 



2 MODEL DESCRIPTION 

The basic mechanisms describing chemical enrichment in 
galaxies can be reduced to infall of gas, metal-rich outflows 
triggered by supernovae, and a star formation prescription. 
Ferreras & Silk (2000, hereafter FSOO) described one such 
simple model for an exponentially decaying infall rate of pri- 
mordial gas in terms of five parameters: flow rate, timescale 
and delay of infall, ejected fraction in outflows, and star 
formation efficiency. Monotonically decreasing star forma- 
tion rates (SFRs) always encounter the so-called G-dwarf 
problem. In FSOO the problem was avoided by assuming 
the metallicity effect was predominant with regard to the 
spectrophotometric output, thereby enabling us to convolve 
simple stellar populations with different metallicities given 
by the chemical enrichment equations but with an average 
common age. 

The G-dwarf problem arises when comparing the pre- 
dictions of closed box models with the observed paucity 
of low-metallicity low-mass stars in the halo of our galaxy 
(Van den Bergh 1962). However, this problem extends to all 
galaxies and morphologies (Worthey et al. 1996) since any 
model assuming a monotonically decreasing infall rate over- 
produces low-metallicity, low-mass stars. These models can 
be readily ruled out for bright cluster early-type galaxies be- 
cause the convolution in age and metallicity of simple stellar 
populations with the chemical enrichment tracks predicted 
by these models yield U — V or V — K colours which are 
significantly bluer than the observed values, U — V ^ 1.6, 
V-K 3.3 for Coma or Virgo (Bower, Lucey & Ellis 1992). 
The solution to this problem involves either pre-enrichment, 
i.e. shifting the zero-point of the metal abundance by as- 



suming non-primordial infall (e.g. Sansom & Proctor 1998), 
or imposing an initial stage of moderate star formation that 
enriches the interstellar medium, followed by a second stage 
where most of the stars (with non-primordial abundances) 
are formed. Several authors (Elbaz, Arnaud & Vangioni- 
Flam 1995; Chiappini, Matteucci & Gratton 1997) have sug- 
gested two-stage processes to solve this problem, as well as 
a top-heavy initial mass function for the first star forming 
stage that overproduces high mass stars for a prompt en- 
richment of the interstellar medium, motivated by the fact 
that starbursts can be best explained by top-heavy mass 
functions (Chariot et al. 1993). 

In this paper we present a model which avoids these 
assumptions and solves the G-dwarf problem by consider- 
ing infall of primordial gas with a rate given by a gaus- 
sian distribution. This infall combined with a linear Schmidt 
law results in a gaussian profile for the star formation rate 
which is a suitable way of parametrizing a strongly peaked 
starburst. Furthermore, the extension of this profile to disk 
galaxies, and in particular to our own galaxy, accounts for 
the observed metallicity distribution of stars (Rocha-Pinto 
& Maciel 1996, Wyse & Gilmore 1995): the first stars are 
formed in a stage with a low SFR, but their contribution to 
the enrichment of the IGM, along with an increasing SFR 
with time yields a large population of stars with metallici- 
ties close to solar. Once the IGM has reached a very large 
metallicity — after the peak in the SFR — a low forma- 
tion rate keeps the tally of stars with a high metal content 
as low as observed. In this scenario, the difference between 
early-type and late-type galaxies would amount to an extra 
component in the SFR for late-type systems, which continue 
to form stars at a slow rate (Ferreras & Silk, in preparation). 
Hence, the work in this paper uses the same chemical enrich- 
ment equations (described in FSOO) based on the formalism 
described by Tinsley (1980), with just four parameters de- 
scribing the star formation process, namely: 

• Star Formation Efficiency (Coff): The star formation 
rate is modelled by a linear Schmidt law, whose propor- 
tionality constant determines the formation efficiency. This 
parameter is a phenomenological approximation to the com- 
plex physics underlying star formation. Its inverse represents 
the timescale in Gyr over which the SFR is extended for a 
sharply peaked infall of gas. In the Instantaneous Recycling 
Approximation (IRA), the SFR {tp{t)) and the infall rate 
{f{t)) are related by: 

/•oo 

Ht) = CoS / ds/(t-s)e-^''"', (1) 
Jo 

" Ceff [1- (l-Bout)i?]' 

Bout is the ejected fraction of gas in outflows — defined 
below — and R is the returned fraction. Low values of Ceff 
extend the process of star formation, generating a larger age 
spread of the stellar populations. 

• Gas Outflows (-Bout): Part of the metal-enriched gas 
in the IGM is heated by supernovae and ejected from the 
galaxy, decreasing the effective yield. This will be modulated 
by the gravitational potential well of the galaxy, so that 
bright and massive galaxies have lower ejected fractions and 
so, higher metal abundances (Larson 1974, Arimoto & Yoshii 
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1987). -Bout represents the fraction of gas being ejected at 
each timestep. 

• Formation redshift (zf): This parameter describes the 
epoch — t{zp) — at which the infall rate (or roughly the 
SFR) was maximum. 

• Infall timescale (r/): Determines the duration of infall, 
so that the infall rate can be written: 

\2" 



f{t,Tf,ZF) 



exp 



(3) 



For a given set of these four parameters (Ccff,-Bout, 
ZF,Tf) the equations give a chemical enrichment track which 
is used to convolve the simple stellar populations from the 
models of Bruzual & Chariot (2000) both in age and metal- 
licity, to obtain a spectral energy distribution for a given 
redshift. Motivated by recent results of the angular power 
spectrum of the Cosmic Microwave Background (Melchiorri 
et al. 1999) , we assume a flat cosmology with a cosmological 
constant {Qm = 0.3, = 0.7, Ho = 60 km s"^Mpc"^), 
although figure 1 shows the result for an open cosmology 
(JIa = 0) in the range of infall parameters (r/ and zp)- The 
fiducial Initial Mass Function (IMF) chosen is a hybrid one 
between the Scalo (1986) and the Salpeter (1955) IMF with 
mass cutoffs 0.1 < M/Mq < 60. We have adopted the be- 
haviour of the Scalo IMF for low stellar masses (M < 2Mq) 
and a Salpeter IMF for the high mass end. The shallow 
slope of the former for low masses better explains the ob- 
servations, whereas the steeper slope at high masses of the 
Salpeter IMF seem to explain better the stellar populations 
in starbursts. For comparison purposes we show the differ- 
ence between a Salpeter IMF and our mass function when 
calculating the correlation between mass-to-light ratio and 
mass (§5). 

Broadband colours are computed from the spectral en- 
ergy distribution obtained for a given star formation his- 
tory described by these four parameters (Cefi, Bout, zf,t/). 
The colour-magnitude relation observed in local clusters is 
used as a constraint, so that a one-to-one mapping is as- 
sumed between colour and absolute luminosity. We have 
used the observations of Coma early-type galaxies (ellipti- 
cals and lenticulars) by Bower et al. (1992), whose colour- 
magnitude relation has a scatter (±0.05 mag) as small as 
the uncertainties expected from population synthesis mod- 
els (Chariot, Worthey & Bressan 1996). Hence, for a given 
U — V colour obtained from the spectral energy distribution 
we infer an absolute luminosity: 



Mv 



{U-V)+ 0.3830 



(4) 



0.0871 

This correlation was obtained by a robust two-stage linear fit 
technique which computes a first estimate of slope and zero 
point using absolute deviations (e.g. Press et al. 1992) and 
then applies a least-squares fit to the resulting sample after 
culling points that deviate more than a given threshold. This 
method prevents outliers from contributing significantly to 
the final slope and zero point. In this case the number of 
outliers was 4 out of 36 galaxies. Once the V-band absolute 
luminosity is found, the total stellar mass is obtained using 
the predicted mass-to-light ratio. 
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Figure 1. Allowed values of infall timescale {tj) and formation 
redshift (zp) taking into account colour constraints for the bright- 
est galaxies in local and distant clusters: The z = 0.02 constraint 
requires U — V and y — if to be correlated with a difference of 
just ±0.05 mag as observed in the Coma cluster (Bower et al. 
1992); the z = 0.55 contour requires the F555W— F814W colour 
to be redder than 2.4 (C10016±16, Ellis et al. 1997), the z = 1.26 
contour imposes a colour R — K > 5.7 as observed by Rosati et al. 
(1999) in cluster RX J0848. 9+4452. The cosmology used assumes 
a matter density = 0.3 and Hg = 60 km Mpc The 
left and right panels show the contours for an open (S^a = 0) and 
flat (Qa = 0.7) cosmology, respectively. 



3 CONSTRAINTS FROM BRIGHT CLUSTER 
GALAXIES 

The colour of the brightest galaxies in clusters at several red- 
shifts is one of the strongest and most robust constraint that 
can be imposed on the model parameters. Rest frame U — V 
colour is very sensitive to age (as well as metallicity) . The red 
U — V colours of the brightest early- type systems observed in 
Coma {U — V ~ 1.6) need a dominant population of old stars 
with a metallicity around solar or higher. Hence, large infall 
timescales (r/) or low formation redshifts are readily ruled 
out since they predict bluer U—V colours. Furthermore, 
observations of clusters at moderate redshifts (e.g. Ellis et 
al. 1997; Stanford, Eisenhardt & Dickinson 1998; Dressier 
et al. 1999) have found a colour-magnitude relation consis- 
tent with passive evolution, imposing again an old age and 
high metal content. This motivates the need to use a low 
ejected fraction (Bout ~ 0) as well as a high star formation 
efficiency (Coff <; 5) in order to reproduce the spectropho- 
tometric properties of these galaxies. Figure 1 shows the 
constraint on infall parameters for the brightest galaxies us- 
ing a grid of models with Bout = and Ccft = 10. The 
~ constraint imposes a colour— colour correlation be- 
tween U—V and V — K as observed in Coma and Virgo 
(Bower et al. 1992) with a scatter of ±0.05 mag. Notice that 
as long as the bursting stage is short enough (r/ ^ 1 Cyr), 
formation redshifts as low as ~ 1 can be accomodated 
with the photometric data (Ferreras et al. 1999). In fact, a 
simple stellar population (which is the extreme case of the 
model described in this paper as r/ 0) gives a z = color 
of U — V 1.5 for a stellar population with solar metal- 
licity and formation redshift zf = 1. This corresponds to 
a bright {Mv ~ —21.6) cluster galaxy. The z ~ 0.5 con- 
straint imposes a colour F555W— F814W (which maps into 
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restframe U — V) redder than 2.4 for the brightest cluster 
galaxies, as observed in C10016+16 (2 = 0.55, Ellis et al. 
1997 ). The high redshift constraint requires R — K ^ 5.7 
as observed in cluster RX J0848.9+4452 {z = 1.26, Rosati 
et al. 1999). These constraints rule out models with infall 
timescales t/ ^ 1.5 Gyr, which reduces the main star for- 
mation episode to less than ~ 3 Gyr. The formation redshift 
is less restricted, allowing the possibility of zf ^ 3, or even 
zf ^ 2 with a flat cosmology with a cosmological constant 
(!^A = 0.7). 

The constraint on the allowed volume of parameter 
space with respect to infall and its subsequent process of 
star formation shown in figure 1 is rather robust as it takes 
into account an observable with a very small uncertainty, 
namely the broadband colours of the most luminous galax- 
ies. Any significant process of star formation at later stages 
will readily translate into bluer colours which are ruled out 
by observations. The modelling of star formation rate as a 
gaussian profile allows the extension of single burst mod- 
els to processes where several bursts of star formation take 
place, as in a "gaseous" merger scenario. The spectrophoto- 
metric predictions of a multi-burst scenario can be mimicked 
by a single gaussian infall rate with a longer time scale (Fer- 
reras & Silk 2000b). Hence, if we consider RX J0848. 9-1-4452 
as a representative cluster which will evolve into something 
like Coma or Virgo, then any hierarchical merging scenario 
should yield gaseous mergers only at epochs earlier than 
zf ~ 3. Merging stages at lower redshifts should not trigger 
any significant star formation, for example involving gas- 
poor progenitors and hot g hinted at by observations 
of red merging galaxies in cluster MS1054 — 03 (z = 0.83 
Van Dokkum et al. 1999). 

It is worth emphasizing at this point that the bright- 
est galaxies we are dealing with in this model represent the 
brightest "non-peculiar" systems. Hence, the very brightest 
cluster galaxies (BCGs) such as cDs should be excluded. 
In a recent paper, Aragon-Salamanca, Baugh & Kauffmann 
(1998) analysed the evolution with redshift of a sample of 
BCGs in 25 clusters in the redshift range < 2: < 1. They 
found no evolution, or even a fading of the absolute K-hand 
luminosity which can be explained by a significant change (a 
factor of 2) in the stellar mass of these galaxies between red- 
shifts z — 1 and 0. In the formalism of our model, the reddest 
— thereby brightest — galaxies are obtained in the corner 
of parameter space associated with high star formation effi- 
ciency (Ceff ~ 5 — 10). This implies the star formation rate 
will be approximately given by the infall rate 4'{t) ~ /(i). 
For the allowed values of infall parameters (t/, zf) shown in 
figure 1, the change in stellar mass between redshifts 2 = 1 
and is very small, and so this model would not be compati- 
ble with BCGs. This type of galaxy requires a strong merger 
rate at late times (2 ~ 1). Furthermore, the mergers should 
undergo no star formation, involving progenitors with either 
very little gas or hot gas. On the other hand, the evolution 
of the bulk of cluster early-type systems was explored in a 
sample of 38 clusters in a similar redshift range (0.1 < 2 < 1) 
by analysing the near-infrared luminosity function (De Pro- 
pris et al. 1999). A significant positive luminosity evolution 
was found for the Schechter parameter Kt, which traces the 
luminosity of the brightest galaxies. The evolution of is 
found to be compatible with the passive evolution of a sim- 
ple stellar population, which is consistent with the models 




Figure 2. Mass contours — labelled in terms of log M/Mq in 
a 2D parameter space spanned by -Bout and Cgff, according to 
a model with zp = 5 and Tf = 1 Gyr (solid lines) or Tf = 0.5 
Gyr (dotted lines). The only region of parameter space allowed for 
the most massive galaxies requires a high star formation efficiency 
(Ccfl) as well as a low ejected fraction (-Bout). Also shown in the 
figure are the trajectories of a "-Bout sequence" , a "Ceg sequence" 
as well as a mixed "-Bout + Ccff sequence" as defined in the text. 
Notice the mass contours do not differ much between models with 
infall timescales of Tf = 0.5 and 1 Gyr for high SF efficiencies (top 
part of the figure), as expected, since high values of C'cS minimize 
the age spread of the stellar populations. 

presented in this paper for the range of infall parameters 
shown in figure 1. 



4 EXPLORING PARAMETER SPACE 

Once we have fixed a region in parameter space for the 
brightest "non-peculiar" cluster galaxies, we can extend the 
model to all early-type systems by exploring the parameters 
describing star formation efficiency and fraction of ejected 
gas in outfiows. Figure 2 shows a contour of stellar masses 
for two realizations of infall {zp — 5, Tf — 0.5 Gyr, dotted 
lines) and {zp = 5, r/ = 1 Gyr, solid lines). The only con- 
straint imposed is the colour-magnitude relation in Coma 
in order to relate colour — obtained from the star forma- 
tion history described by a chosen set of parameters (-Bout, 
Ccft, Tf, Zf) — and luminosity, or stellar mass. The contours 
change slightly with infall parameters but the qualitative be- 
haviour is unchanged. We need bot-h a high star formation 
efficiency and a low ejected fraction in outfiows to account 
for the most massive galaxies. Lowering the SF efficiency will 
generate a significant spread in stellar ages that will yield 
colours that are too blue compared with the observations. 
On the other hand, a higher ejected fraction will decrease 
the effective yield, so that the average metallicity will be 
lower, blueing the predicted colours with respect to the ob- 
servations. Hence, regardless of the choice of parameters, 
we find the brightest cluster galaxies to lie in a well-defined 
corner of the 4-dimensional parameter space explored in this 
paper. Lower mass galaxies are harder to identify because 
of a degeneracy between outfiow rates and star formation 
efficiency. The blueing of galaxies as we progress down in 
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Figure 3. Distribution of stellar ages and metallicities according 
to the three sequences defined in the text. The top panels show 
the average age (ieft) , and metallicity (right) , whereas the bottom 
panels show the standard deviation of the distributions. Notice a 
"Sout sequence" generates a pure metallicity sequence for cluster 
ellipticals, with all galaxies being roughly coeval, in agreement 
with the observations in the Fornax cluster (Kuntscher & Davies 
1998) The closest model to a pure age sequence is a "Cog se- 
quence", where a range of star formation efficiencies result in a 
significant average age spread with respect to luminosity. How- 
ever, there is also an important range of metallicities. A pure age 
sequence is ruled out because of chemical enrichment. 



luminosity along the colour-magnitude relation can be ex- 
plained either by decreasing the average metallicity of the 
stellar populations (i.e. a range of ejected fractions, Bout) or 
by increasing the fraction of younger stars (i.e. a range of SF 
efficiencies, Ces)- The real picture will likely involve a range 
both in Bout and Ceff. At low redshift it is hard to disen- 
tangle the contribution from these two parameters, mainly 
caused by the age- metallicity degeneracy (Worthey 1994). 
However, these two alternative mechanisms will predict dif- 
ferent evolutions with lookback time: as we go to higher 
redshift, the sequence driven by SF efficiency (hereafter C'ctt 
sequence) will show stronger evolutionary scars because of 
the larger scatter in age, whereas the sequence driven by 
outflows (hereafter Bout sequence) will be unaffected until 
the observations go back to the epoch of star formation. 

Figure 3 shows the predicted distributions of stellar ages 
(left) and metallicities (right). The average (top) and stan- 
dard deviation {bottom) of these distributions are shown for 
three sequences: Bout,Ccfr and Bout+CeS driven by outflows, 
star formation efficiency and a mixture of both, respectively. 
All of these models have r/ = 1 Gyr, zf — 5; the Bout se- 
quence fixes Ccff = 5; the Cotf sequence flxes Bout ~ 0; 
and the mixed Bout + Ceff sequence imposes a correlation 
between these two parameters. We expect the star forma- 
tion efficiency to be higher in more massive galaxies (where 



the ejected fraction should be lower), hence we have chosen: 
log Ccff — 1 — 3 -Bout, as an illustrative example. An outflow- 
driven sequence generates a pure metallicity sequence where 
the stellar populations of all galaxies are coeval. This is the 
model favored by the analysis of the colour-magnitude rela- 
tion by Kodama et al. (1998) for a sample of 17 clusters in a 
wide redshift range (0.3 < z < 1.3). In a Bout sequence the 
colour-magnitude relation is due to a pure mass-metallicity 
correlation, so that the brightest galaxies have higher metal 
abundances. On the other hand, a Ccff or a Bout +Ccft model 
represents mixed age-|-metallicity sequences as it generates 
a significant age spread of the stellar populations as well as 
a metallicity range. A pure age sequence explored in other 
papers (e.g. Kodama & Arimoto 1997) is ruled out when 
including chemical enrichment for any reasonable assump- 
tion about the star formation history of the galaxy. The low 
metallicity expected in the fainter galaxies can therefore be 
explained either by a lower effective yield caused by gas out- 
flows, or by a low star formation efficiency that slows down 
the process of enrichment. Nagashima & Gouda (1999) men- 
tion the effect of the UV background radiation as a possible 
alternative to supernovae winds to suppress chemical enrich- 
ment. However, in this scenario it seems hard to motivate the 
correlation found with respect to luminosity (e.g. a stronger 
feedback should be expected from the UV background for 
the fainter galaxies). Figure 3 shows that a Ccff sequence 
implies a ~ 3 Gyr age difference between the brightest sys- 
tems and Mv ~ — 16 galaxies with a larger age spread for the 
fainter galaxies, as expected for a system with a lower star 
formation efficiency. This model would thus be consistent 
with the smaller age and metallicity scatter found in Coma 
ellipticals (j0rgensen 1999). This age difference is negligible 
for local clusters, since the predicted youngest stars are still 
too old (~ 8 Gyr) to be detected spectrophotometrically. 
However, predictions between a coeval sequence and a Ccff 
sequence at lookback times ~ 8 Gyr (i.e. z ^ 1) differ signifi- 
cantly as shown below (§5). A Bout sequence predicts a small 
age and metallicity spread in faint ellipticals. In this case the 
observed photometric scaling relations are just driven by the 
average metallicity. 



5 MASS-TO-LIGHT RATIOS: STELLAR 
VERSUS OBSERVED 

The evolution of the stellar mass-to-light ratio is a suitable 
probe for the analysis of the ages of the stellar populations 
in galaxies because of its weak dependence on metallicity 
(especially when measured in NIR passbands). However, ob- 
served M/L ratios carry large uncertainties and systematic 
offsets as they are computed from central velocity disper- 
sions (cro), sizes (re) and surface brightnesses (Ee), requir- 
ing a prescription for the structure in order to infer masses 
and luminosities. From a simple dimensional analysis, we 
can write the mass and luminosity of a galaxy as: 



M = aa^re 
L ■ 



(5) 
(6) 



The proportionality constant {3 only depends on the way the 
surface brightness is defined with respect to the effective size. 
However, a depends on the structure of the galaxy, which re- 
quires a dynamical model (e.g. King 1966; Bender, Burstein 
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Figure 4. Model predictions of the correlation between stellar 
mass-to-light ratio and mass in two different passbands: V band 
(top) and K band {bottom). Notice the NIR observable is much 
more sensitive to age and thus changes significantly with infall 
timescale. The tj = 1 Gyr case also displays the predictions for 
a Salpeter IMF, rather high in comparison with observations. 
The triangles show Coma cluster ellipticals from Mobasher et 
al. (1999), whereas the dots show ellipticals in a sample of 11 
clusters extracted from Pahre (1999). The predictions — which 
just give stellar masses and mass-to-light ratios — have been 
transformed to total mass and M/L using a non-linear corre- 
lation between total mass (inferred by observations), and stellar 
mass (given by population synthesis models): Af-poT °c JV/g^^. 
The grid comprises a set of Bout sequences (thin) with fixed 
Ceff = {0.1,0.5, 1,5, 10} and a set of C^ff sequences (thick) with 
fixed Bout = {0,0.2,0.4,0.6,0.8}. 

& Faber 1992), usually assumed to be invariant with regard 
to galaxy mass or size (i.e. homologous). Furthermore, a 
comparison between observations and model predictions de- 
pends on a prescription that relates stellar and total matter. 
Despite all these caveats, mass-to-light ratios are one of the 
best candidates for breaking the age-metallicity degeneracy. 

Figure 4 shows the predicted M/L ratios in two pass- 
bands: optical [V) and NIR {K) for a grid of Bout (thin) and 
Ceft (bold) sequences. A Bout + Cotf mixed sequence would 
be represented by a trajectory between the grid spanned 
by these lines. The results are shown for three different in- 
fall timescales r/ = 0.5, 1, 1.5 Gyr. The solid triangles are 
the observations of Coma cluster galaxies from Mobasher 
et al. (1999), whereas the dots represent ellipticals in 11 
clusters from the sample of Pahre (1999). The slope of the 
correlation between mass-to-light ratios and masses is 0.24 
in the V band and 0.14 in K band, whereas purely stellar 
M/L ratios yield slopes around 0.1 and 0.0 in the V and K 
bands, respectively. This mismatch cannot be related to a 
different stellar population. Being a strongly age-dependent 
observable, the only way to generate the observed slopes 
in M/L would require a stellar age spread with respect to 



galaxy mass whose restframe U — V colour would signifi- 
cantly disagree with the colours in moderate redshift cluster 
ellipticals. As a simple check to confirm this, we used simple 
stellar populations (SSP, i.e. no age or metallicity spread). If 
we relate the brightest galaxies {My ~ —22.5) to a SSP with 
an age of 12 Gyr and a super-solar metallicity Z = 1.5Zq, 
the required age for the SSP for a Mv = — 19 galaxy which 
would explain the observed slope M/Lv oc M"'^** should be 
f = 4 Gyr (Z = Zq) or t = ^ Gyr {Z = Zq/2). These 
two possibilities yield U ~ V colours of 1.20 and 1.07 re- 
spectively, so that the latter is barely consistent with the 
observed broadband photometry in nearby clusters (Bower 
et al. 1992) . However, this hypothesis fails when considering 
moderate redshift clusters: for instance, cluster C100164-16 
{z ~ 0.5) has a well-defined red-envelope in a large range of 
luminosities (Ellis et al. 1997), yet the lookback time (~ 6.2 
Gyr in our adopted cosmology) would be larger than the 
expected age for the stars in the fainter galaxies ! 

Hence, the only plausible way to solve this slope mis- 
match is by imposing a non-linear correlation between total 
mass (including dark matter) and stellar mass. A simple 
power law M oc (FSOO) has been used in figure 4 to 

transform predicted stellar masses and mass-to-light ratios 
to total masses and M/L. The agreement is good in both 
bands, which means the observed slope difference in the fun- 
damental plane between passbands can be explained by the 
stellar populations alone. This result agrees with the anal- 
ysis of j0rgensen (1999) who refers to a higher fraction of 
dark matter in the most massive ellipticals in order to ac- 
count for the mismatch between age estimates using M/L 
ratios (which involves the total mass of the galaxies) or HPc 
indices (involving just the stellar component). However, one 
should bear in mind that there are alternative explanations 
for this slope mismatch, such as a systematic variation of 
the initial mass function or the breaking of the homology for 
the structural parameters of early-type systems (j0rgensen 
1999; Pahre et al. 1998). Graham & CoUess (1997) examined 
the effect of a non-homologous light profile in 26 early-type 
systems in the Virgo cluster and found no differences in the 
fundamental plane from estimates using a universal de Vau- 
couleurs profile. 

Notice the scatter of the correlation is not the minimum 
one can obtain when computing the fundamental plane. The 
dependence of the fundamental plane on the observables (ve- 
locity dispersion, surface brightness and galaxy size) is not 
precisely the one inferred for the correlation between mass- 
to-light ratio and mass. This scatter has an observational 
component as well as an intrinsic contribution. The intrin- 
sic part can be related to differences in either the stellar ages 
or the initial mass function of galaxies with the same mass. 
The former can be caused by slightly different formation 
redshifts or infall timescales. Notice the predictions of M/L 
ratios for three different infall timescales in figure 4. A vari- 
ation of 1 Gyr in infall timescale can account for changes 
in log M/L roughly of order 0.1 dex. Furthermore, a non- 
universal IMF will add a variation around 0.2 — 0.3 dex. 
The top-center panel displays the predictions for our fidu- 
cial "hybrid" IMF and a Salpeter IMF with the same mass 
cutoffs (0.1 < M/Mq < 60). a different cutoff or slope at 
the high mass end will not change the result very much since 
these changes only have an effect on metallicity. However, 
changes in the low mass end of the IMF will significantly al- 
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and zero point {bottom panels) of the correlation between 
mass-to-light ratio and mass for a range of infall timescales 
(left) or formation redshifts (right). The points show the 
observations of a few clusters in restframe B band: Coma 
{z = 0.02), CI 1358-^62 {z = 0.33), CI 0024-f 16 (z = 0.39), 
MS 2053-h03 (z = 0.58) and MS1054-03 (2 = 0.83) from 
the compilation in Van Dokkum et al. (1998). Any sequence 
(-Bout, Ccft, -Bout -I- Ccff) will predict a decrease of the zero 
point as the stellar populations get younger at higher red- 
shift, although the zero point decreases faster for a Ccft se- 
quence, for which the average stellar population is younger. 
Current observations of moderate redshift cluster ellipticals 
allow an estimate of the zero point but not of the slope with 
enough accuracy. For instance, the study of j0rgensen et al. 
(1999) with clusters at low and moderate redshifts seem to 
indicate a steepening of the correlation between M/L and 
mass, thereby favouring a Ccff sequence. However, selection 
effects could actually mimic such behaviour. Forthcoming 
spectroscopic observations of cluster ellipticals will enable 
us to estimate this slope, thereby allowing us to determine 
the importance of outflows and star formation efficiency in 
the evolution of cluster early-type galaxies. 



Figure 5. Predicted evolution with observed redshift of the slope 
(top) and zero point (bottom) of the correlation between mass- to- 
light ratio and mass in the B band, in models driven by star for- 
mation efficiency (Cgjf , thick lines) or outflows (-Bout, thin lines). 
rjg is defined as the difference between the slope in log M/Lg 
vs log -M at 2: > and z = 0. The data points in the bottom 
panels show observed shifts in the zero point for clusters: Coma 
(2 = 0.02), CI 1358-1-62 (z = 0.33), CI 0024-1-16 (z = 0.39), MS 
2053-1-03 (z = 0.58) and MS1054-03 (z = 0.83) from the compi- 
lation in Van Dokkum et al.(1998). 



ter the mass-to-light ratio. Our "hybrid IMF" behaves like a 
Scalo function at low masses, whereas a Salpeter mass func- 
tion produces a large fraction of low mass stars, thereby 
increasing the stellar M/L ratio. 

The different age distributions obtained for the -Bout 
and Ccff sequences generate different predictions for the evo- 
lution of the slope of the correlation between stellar M/L 
ratio and mass. We have shown above that a -Bout sequence 
(driven by outffows) is equivalent to a metallicity sequence. 
Hence, a very small slope change is expected since mass- 
to-light ratios are not very sensitive to metal abundance. 
Furthermore, the decrease in M/L will be greater at higher 
metalhcities (i.e. for the brightest galaxies), so that the cor- 
relation will get ffatter. In order to quantify the slope change, 
an index rjx is defined as the slope change between redshift 
2; > and z — using mass-to-light ratios in the X pass- 
band, namely: 



Vx(z) 



A log M /Lj. 



A log M 



A log M/L^ 



A log A'/ 



(7) 



Hence, a -Bout sequence predicts small negative values for 
rjx{z)- On the other hand, a Ccff sequence is driven by age, 
so that the faintest galaxies — which will have a younger 
stellar population — will decrease their M/L ratios with 
redshift faster than the brighter (and older) galaxies. This 
corresponds to a more significant and positive value for 
rix(z). Figure 5 shows the evolution of the slope (top panels) 



6 DISCUSSION 

A simple phenomenological treatment is described in this 
paper, where the mechanisms underlying the evolution of 
galaxies are reduced to a set of four parameters. The star 
formation rate is assumed to follow a linear Schmidt law 
whose proportionality constant is used to describe a varying 
star formation efficiency (Ccff). The supply of primordial 
gas fuelling star formation is controlled by gaussian infall 
characterised by the epoch at which the infall rate is max- 
imum (zf), and the width of the gaussian profile gives a 
characteristic infall timescale (r/). The model is allowed to 
eject a fraction (-Bout) of the enriched gas, thereby lower- 
ing the effective yield. A first stage in this analysis involves 
finding a suitable pair of infall parameters (zF,Tf) which 
are capable of generating the restframe U — V colours of 
the brightest cluster galaxies. The colour constraints im- 
posed by the reddest (and brightest) galaxies in clusters at 
moderate and high redshift allow us to discard long infall 
timescales and recent star formation epochs. One could ar- 
gue that the constraint on the star formation history of the 
brightest systems need not be the same as for less massive 
ellipticals. However, the current most plausible scenarios for 
galaxy formation assume either a simultaneous process of 
star formation regardless of galaxy mass, or a hierarchical 
structure where the most massive galaxies might have un- 
dergone the latest bursts of star formation. Furthermore, the 
possibility of an "inverted-hierarchical" scenario can still be 
accomodated in this model, as long as the constraints im- 
posed (zf <; 3,r/ <; 1 Gyr) are held even for the less massive 
galaxies. 

Out of the four parameters considered in the model, 
we have found that — within the framework of this model 
— only the SF efficiency (Ccff ) and the ejected fraction in 
outflows (-Bout) help determine the mass sequence in cluster 
early-type galaxies. Infall parameters (zF,Tf) do not appre- 
ciably change much the output unless very recent stages 
of star formation are included. However, this will result in 
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restframe U — V colours that are in contradiction with ob- 
servations. The efficiency parameter generates a significant 
spread in the age distribution of stars, ahhough mixed with 
a metalhcity range (any model with a reasonable IMF must 
include this range of abundances). Alternatively, a range 
of outflow rates result in a range of metallicities with no 
significant spread in ages. Both mechanisms are degener- 
ate in local clusters because the average age of the stellar 
populations predicted for any model at z — are too old 
to be able to disentangle the effects of age and metallicity. 
However, the predictions of age-sensitive observables at high 
redshift differ noticeably for these two sequences. Unfortu- 
nately, present data is not capable of ruling out one model 
against the other, but the continuing flow of data from clus- 
ters at moderate-to-high redshift will eventually enable us to 
break this degeneracy. A sequence driven by efficiency (Ccff 
sequence) predicts a steepening of the slope of the funda- 
mental plane, whereas a -Bout sequence — driven by outflows 
— predicts no change or a slight decrease of this slope. High 
precision observations of the dynamical and spectrophoto- 
metric properties of cluster galaxies at high redshift will con- 
firm the importance of either age or metallicity in the mass 
range of cluster ellipticals. 
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